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ABSTRACT 

We present NIR observations of a sample of mainly interacting ULIRGs, comprising 
H and K band spectra. Our main aims are to investigate the power source of these 
extremely luminous objects and the various excitation mechanisms of the strong molec- 
ular hydrogen emission often seen in such objects. Broadened emission lines were only 
detected in one object, IRAS 23498, consistent with previous results for this galaxy. 
The [SiVI] emission line was detected in IRAS 17179 and IRAS 20210, both classified 
as Sy2s. Two of the sample were unclassified, IRAS 00150 and IRAS 23420, which ex- 
hibit neither [SiVI] emission or broadened HI emission. However this does not rule out 
the presence of an obscured AGN. Analysis of the molecular hydrogen emission showed 
that the major excitation mechanism for most objects was thermal. Modelling of the 
more luminous objects indicate for IRAS 20210 10 per cent, and for both IRAS 23365 
and IRAS 23420, 30 per cent of their 1-0S(1) line emission has a non-thermal origin. 

Key words: galaxies: active - galaxies: interactions - galaxies: nuclei - galaxies: 
Seyfert - galaxies: starburst - infrared: galaxies. 



1 INTRODUCTION 

The existence of ultra-luminous infrared galaxies (ULIRGs) 
was highlighted by the Infrared Astronomical Satellite 



many cases the identity of the central source is still uncer- 
tain. 

A large number of ULIRGs have now been i dentified 
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Many have been well studied via optical imaging ( Leech ct 
al. 1994 IClements et al. 1996a| ; [Auriere et al. 19961 iMurphy 
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— 0. Nate that, unless otherwise specified, we have adopted 

H = 50 km s _1 Mpc -1 and q = throughout. In order to 
produce energy output of this magnitude, ULIRGs require 
a substantial energy source. The two main candidates are a 
burst of strong star formation in a central starburst or an 
active galactic nucleus (AGN) with accretion onto a central 
black hole. ULIR Gs are generally very gas and dust rich e.g. 
Rigopoulou et al.( 1996). Unfortunately this gas and dust ob- 



scures the centre of the galaxies in which they reside, ham- 
pering identification of the central source. Direct evidence 
for the presence of AGN is found for some ULIRGs in the 
form of Seyfert-like emission line ratios, high excitation lines 



and sometimes broad permitted-line widths (Veilleux et al 
1995). It has been proposed that the fraction of ULIRGs 



contai ning A GN rises with infrared luminosity (e.g. Sanders 
et al. ( 1988 )). Often evidence for the pr esence of strong star 
formati pn io found (o.g. Gonad ot al. ( 1008 )). Howovor in 



et al. 1996 ;|Suracc ct al. 2000) and spectros copy (Lutz et al 
199S| ; [Kim et al. 1998| ; |Vcillcux et al. 1999^ ). However, data 
of longer wavelengths, such as the near infrared (NIR), have 
the potential to probe deeper into the dust which enshrouds 
these galaxies, since the extinction is much less i.e. Ay ~ 10 
A k . Relative to t he opt ical, less NIR data exists for ULIRGs. 
Goldader et al. (1995) have obtained K band spectra for 
a sample of 13 ULIRGs and 24 luminous infrared galaxies 
(LIRGs). Comparison of the properties of the two sets of 
galaxies highlights a suppression in the dereddened lumi- 
nosity of the Br7 emission, relative to the far-infrared lumi- 
nosity for th e ULIRGs, as compared to the LIRGs. Murphy 
et al. (1996) have studied a sample consisting of 46 galaxies 
from the IRAS 2 Jy survey. They have imaged these galaxies 
in the optical and the NIR, and find a high percentage to 
be interacting (95 per cent). They have also obtained NIR 
spectral data for 33 of t hese objects for wh i ch results has 
recen tly been published ( Murphy ct al. 1999| ; Murphy et al 
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200C| ). Finally, Veilleux et al. ( |1997| ; |1999b| ) have conducted 
a campaign aimed at uncovering hidden broad line regions 
in ULIRGs. They have observ ed in total 6 4 ULIRGs, taken 
from the IRAS 1 Jy sample (|Kim 1995|; [Kim et al. 199&j; 
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Kim & Sanders 1998). Broad Paa is detected for the first 



time in four sources and possibly more, some uncertainty is 
due to the potential contamination from the presence of a 
nearby H2 line. 

We have selected a subs et of galaxies mainly from the 
sample of Leech et al. (1994). These display a cross-section 
of interaction characteristics i.e. ranging from some that are 
obviously interacting to those with no companion within a 
significant distance. Two other ULIRGs were added which 
show evidence of interaction in their morphology. Our final 
sample shows a bias towards interacting galaxies. 

H and K band spectra have been obtained for this sam- 
ple. Using this data we have investigated the central sources 
of the ULIRGs, looking for possible broadened emission lines 
and/or high excitation lines indicative of the presence of an 
AGN. We also study the excitation mechanisms of the lines 
present. This will allow us to determine which processes are 
dominant in particular objects, for example excitation due 
to young stars, supernova remnants (SNR), a central black 
hole or shocks due to interaction. Additionally we investi- 
gate further relationships which have been established previ- 
ously, such as the under-luminous Br7 line issue mentioned 
previously. In addition we have J, H and K band imaging 
f or our sample which is discussed separately in Davies et al. 



( ^OOObj ). 



This paper is organised as follows. First, the sample is 
discussed in § |i| The spectral data are then presented and 
the observations and data reduction are discussed in § || 
The final results are presented and discussed in § ^. The 
summary of our results appears in § ^, and a brief synopsis 
for each galaxy is included in Appendix 



2 THE SAMPLE 

The galaxies in our sample were selected based on their 
extreme infrared luminosities (Lir > 1O 12 L0, using H = 
50 km s _1 Mpc -1 and q = 0) . We note that when using 
the Sanders and Mirabel ( 1996| ) definition of this classifi- 
cation, which uses the higher value of H , two of our sam- 
ple, IRAS 20210 and IRAS 23420, have values for Lir falling 
slightly below 10 12 L . The redshift range of the sample, z 
= 0.06 - 0.21, was chosen to be well suited to ground-based 
NIR observations. For all the objects Paa is shifted into the 
K band, except IRAS 23498 which was only observed in the 
H band in which Pa/3 falls. The [Fell] line at 1.644 /im and 
a number of H2 lines also fall into the K band window. 

Five of the seven objects in our sample are a small sub- 
set of a l arge sample observed in the optical by Leech et 
al. (1994). Their original sample was complied from four 



flux limited (at Sqq = 0.6) su b-samples of the QDO T red- 
shift survey of IRAS galaxies (Lawrence et al. 1999). These 
consist of (1) all known ULIRGs from the No rth G alac- 
tic Wedge (NGW) survey of Lawerence et al. (Il986|), (2) 



all known ULIRGs from the QDOT survey in a specified 
RA and DEC range, (3) other randomly selected ULIRGs 
from the QDOT survey and (4) a number of high luminos- 
ity (5 x 10 11 < -LeOfim/i© < 10 12 ) galaxies from the same 
survey. 

IRAS 20210 and IRAS 23365 were not part of the Leech 
sample but were taken from the literature as additional, 
bright examples of nearby interacting ULIRGs falling within 



the observing constraints. We note that this sample is not 
complete in any sense. 

Table [I] lists the basic properties of the galaxies pre- 
sented in this paper. Galaxy names are given in column (1), 
positions in columns (2) and (3), interaction class in column 
(4), redshift, z, taken from the NASA/IPAC Extragalactic 
Database in column (5), infrared luminosity in column (6) 
and galaxy classification taken from the literature, if known, 
in column (7). References for classifications are given in Ap- 
pendix [a|. 



3 OBSERVATIONS AND DATA REDUCTION 

The spectral observations were performed using the United 
Kingdom Infrared Telescope (UKIRT) on the nights of 7th 
and 8th August 1999. The long slit grating spectrometer, 
CGS4, was used with H and K band filters. The orientation 
of the slit during observations was chosen to be along the 
major axis of the galaxy or across both nuclei if a second 
nucleus was known to be present. Integration times were 
usually of the order of ~ 20 - 40 minutes, per object, in 
both bands. Sky frames were obtained by nodding along the 
slit, thus allowing continuous observation of the galaxies. 
To properly sample the CGS4 spectrum in the configura- 
tion usually used for these measurements, four individual 
spectra were obtained, with the spectrum shifted by half a 
pixel over 2 pixels. This has the additional benefit of allow- 
ing the removal of the effects of bad pixels in a single row. 
These frames were combined to produce an overall spec- 
trum. Spectra were generally extracted using 5 pixels along 
the slit; outside of these few pixels relatively little (~1 per 
cent) or no signal was detected. The resulting spectra were 
fully sampled at a resolving power of 800 (375 km s _1 ) at 2 
/im, with a slit width and pixel size of 0.61 arcsec projected 
on the sky. 

Basic data reduction was performed using IRAF and 
consisted of the following steps; 

(i) Subtraction of bias frames and division by flat frames 
removing the bias level present in the array and pixel-to- 
pixel variations. 

(ii) Subtraction of sky frames from object frames. This 
removed sky lines and created, effectively, two spectra on 
one image, one being positive, the other negative. 

(iii) Transformation of the images using a matrix pro- 
duced from arc images to remove distortion effects often 
found in long slit images. 

(iv) Addition of positive and negative spectra, producing 
a final spectrum. 

(v) Division by the spectrum of a type F (or earlier) star, 
observed at a similar air mass, correcting the spectra for 
instrumental and atmospheric responses. Multiplication by 
a black body at the effective temperature of the standard 
provided the intrinsic spectral shape of the galaxy. Obser- 
vations of flux standards reduced in the same fashion were 
used to obtain the final flux calibration of the galaxy. 

Early type stars were chosen since they do not show 
metal features seen in later type stars which would have 
masked features in the galaxy spectra or produce false emis- 
sion features after division. The stars used do exhibit some 
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Table 1. Details for ULIRGs observed. 



IRAS Name 


RA (J2000) 
l/J 


DEC (J2000) 
( rj .\ 


Interacting? 


z 

l+J 


log i-LK 

V>) 


Comments 
(I ) 


00150+4937 


00hl7m44.8s 


+49d54ml5s 


Yes 


0.148 


12.64 




16487+5447 


16h49m47.2s 


+54d42m32s 


Yes 


0.104 


12.48 


LINER 1 - 2 


17179+5444 


17hl8m54.2s 


+54d41m47s 


No 


0.147 


12.57 


Sy2 2 


20210+1121 


20h23m25.9s 


+Ild31m31s 


Yes 


0.056 


12.19 


Sy2 3 - 4 


23365+3604 


23h39m01.3s 


+36d21ml0s 


Disturbed 


0.065 


12.45 


LINER 5 


23420+2227 


23h44m32.3s 


+22d44m29s 


No 


0.087 


12.09 




23498+2423 


23h52m25.0s 


+24d40m09s 


Yes 


0.212 


12.82 


Sy2 2 



References: (1) Kim et al. ( [1998; ) (21 V cillcux et al. ( |1999e)) (3) Vader et al. ( |1993[ ) (4) Perez et al. 

( 199C|) (5) Veilleux et al. (|l995|). 



Table 2. Bands and slit position angles used for observations. 



Galaxy 



Band 



PA 



Exp. time 
(minutes) 



(IRAS) 


H 


K 


(east of north) 


H 


K 


00150 


yes 


yes 


172 


24 


40 


16487 


yes+ 


yes 


135 


2d 


10 


17179 


yes 


yes 


57 


40 


56 


20210 


yes 


yes 


166 


24 


40 


23365 


yes 


yes 





40 


10 


23420 


yes 


yes 1 " 


12 


40 


48 


23498 


yes 


no 


108 


24 





f Variable weather conditions during this observation makes flux 
calibration uncertain, see main text for more details. 



atomic hydrogen features (e.g. Br7 in the K band) which 
were removed using Lorentzian fits. 



4 RESULTS AND DISCUSSION 

Table presents details of the observations made for each 
object, listing the specific bands, slit position angles and 
on source exposure times used. The final reduced H and K 
band spectra at observed wavelengths are shown merged in 
Figure |l[ Figure ^ shows the merged spectrum, produced 
using DIPSO, for IRAS 20210 at rest wavelength with the 
prominent lines labelled. The spectra have been scaled to 
the strongest line present, Paa. The only exception to this 
is IRAS 23498 which was not observed in the K band. A plot 
of the atmospheric transmission has been included with the 
plots to indicate areas where noise may dominate the spec- 
trum, effectively obscuring lines present in the affected wave- 
band. These data, produced using the program IRTRANS4, 
were obtained from the UKIRT worldwide web pages. Emis- 
sion from a second nucleus was detected for IRAS 00150 and 
IRAS 20210, these spectra have been included together with 
that of the primary nucleus. In both cases the second nucleus 
has much lower flux level. Note that, due to the slit position 
angles used, the secondary nucleus seen for IRAS 16487 in 
the K' b and (|Clcmcnts ct al. 1996a| ), and IRAS 23498 in the 
optical ( Leech et*aL*T994| ) did not fall within the slit. 

A useful check on the validity of the flux calibration 
was possible since the H and K bands were reduced indepen- 
dently for each object. The agreement between the bands is 
generally good, and lies within 20 per cent. Variability in 
weather conditions made flux calibration uncertain on the 



two occasions indicated in Table |^. In these cases the band 
for which the calibration was uncertain was scaled to match 
the other band, for the same galaxy, for which we are confi- 
dent of the flux calibration. 

Table ^| lists the line fluxes and equivalent widths mea- 
sured using the separate H and K band spectra. In order 
to accurately measure line fluxes a continuum was fitted to 
each spectrum and subtracted. Points free of known spec- 
tral features were used. However, in cases where significant 
variability in atmospheric conditions made removal of atmo- 
spheric effects uncertain, such as for IRAS 23365 in the H 
band, fitting of a continuum was somewhat subjective. In 
order to assess the magnitude of this effect, several low or- 
der fits were tried, subsequent line fluxes varied by less than 
~10 per cent. Due to the overlap in wavelength between 
the H and K bands some lines appear in both spectra e.g. 
Paa is visible in both bands for IRAS 20210, IRAS 23365 
and IRAS 23420. When this is the case the final values for 
flux and equivalent width quoted are averages, except in the 
cases mentioned above where flux calibration is uncertain 
due to variable weather conditions, when the more reliable 
measurement was adopted. 

In the majority of the spectra a blue 'shoulder' is evi- 
dent on the Paa line (A1.876/^m). The origin of this will be 
discussed further in § 4.3. All fluxes and equivalent widths 
quoted in Table [| for Paa include this 'shoulder'. Of note is a 
line blend consisting of contributions from BrS (A1.945/im), 
H 2 1-0S(3) (A1.958Aim) and [SiVI] (A1.963/im), prominent 
in the K band for all objects. To facilitate flux measure- 
ments these lines were deblended using the routine 'splot' 
from IRAF and also using DIPSO. 

Note that Pa/3 was observed in the H band for 
IRAS 23498. Although not given in Table | this line is dis- 
cussed further in 8 i.c . 



4.1 General features of the spectra 

As can be seen the K band spectra are much richer in lines 
than the H band for the redshifts of the galaxies observed. 
In the K band most of the galaxies observed show several 
hydrogen recombination lines, most notably Paa and Br7. A 
number of H2 excitation lines are seen, especially numerous 
for the brighter galaxies i.e. IRAS 20210. Other lines visible 
are Hel at 2.059^m and [Fell] at 1.644/xm. [SiVI] at 1.963/im 
is observed for IRAS 17179 and IRAS 20210. 

Several CO bandheads longward of 2.3/im can be seen in 
two of the galaxies observed at this wavelength: IRAS 20210 
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Figure 1. Merged H and K band spectra of flux at observed wavelength in units of I0~ 12 erg s — 1 cm -2 fim -1 . Each spectrum has been 
scaled to the strongest line present, Pact. This is not the case for IRAS 2349 observed only in the H band. The vertical dashed lines show 
the overlap between the H and K band observations. The uppermost plot shows the atmospheric transmission for the same waveband 
in order to indicate wavelengths where atmospheric effects are important. These data, produced using the program IRTRANS4, were 
obtained from the UKIRT worldwide web pages. For line identifications see Figure |. © 0000 RAS, MNRAS 000, 000-000 
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Figure 1 - continued 



2.2 



2.4 



and IRAS 23365. These are not obvious in Figure [l] due to 
the scaling used. See § |4.8| for further discussion. 

The H band shows few features for the majority of 
galaxies observed, with the exception of IRAS 23498 which 
shows a clearly broadened Pa/3 line with FWHM ~ 5000 km 
s _1 . This is discussed further in S 4.3. 



4.2 Extinction 

Values of internal reddening were calculated using the NIR 
recombination lines, measured here, and optical lines taken 



from the literature, available for four of the seven galaxies 
in this sample. Table ^| lists the NIR - optical line combina- 
tions used and extinction values found for each. The method 



used follows that of Veilleux and Osterbrock (1987), using 
the reddening c urve for the optical and NIR parametrised by 
Cardelli et al. (1989). The intrinsic Hq/H/3 ratio was taken 
to be 2.85 for starburst galaxies (assuming Case B recombi- 
nation) and 3.1 fo r AGN (having enhanced Ha d ue to colli- 
sional excitation (Veilleux and Osterbrock 1987)). Galactic 
extinction for each galaxy is shown in Table ^j. These val- 
ues are small and can be assumed to be negligible, even for 
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Wavelength (fJ-m) 



Figure 2. Merged H and K band spectrum of IRAS 20210 with the main lines observed labelled. The scaling has been chosen so as to 
enhance the less luminous lines seen. Note that atmospheric transmission is particularly poor for the wavelength range 1.7 — 1.85 /an 
and most of the 'features' in this region are noise. 



Table 3. Fluxes and equivalent widths for emission lines seen in spectra. 



Galaxy 


[Fell] 


Br!) 


H 2 


Pan 


H 2 


BrS 


H 2 


[SiVI] 


H 2 


Hel 


H 2 


Br- 


H 2 


H 2 


(IRAS) 


1.644 


1.818 


1.836 


1.876 


1.892 


1.945 


1.958 


1.963 


2.034 


2.059 


2.122 


2.166 


2.223 


2.248 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


00150a 








0.671 




0.016 






O.Oll 31 






0.017 














92.6 




8.4 






6.7 






13.8 






00150b 








0.108 






























50.1 






















16487+ 


0.256 




0.074 


0.938 


0.037 


0.067 


0.160 




0.058 


0.040 


0.182 


0.06 








0.4 




17.0 


216.6 


9.1 


17.7 


42.6 




17.3 


12.6 


62.0 


20.7 






17179 


0.050 




0.060 


0.773 






0.146 


0.110 


0.022 x 


0.018* 


0.092 


0.056 








19.8 




6.9 


90.9 






18.9 


14.3 


2.9 1 


2.4* 


13.0 


7.9 






20210a 


1.237 


0.144 


0.382 


9.371 


0.317 


0.631 


1.241 


0.860 


0.317 


0.404 


0.943 


0.686 


0.236 


0.121 




54.7 


7.1 


18.3 


46.2 


15.9 


32.4 


63.4 


42.8 


16.6 


21.3 


50.2 


36.9 


12.9 


6.8 


20210b 






























23365 


0.163 






2.266 


0.096 


0.068 


0.186 




0.073 


0.068 


0.173 


0.157 


0.051 






7.4 






115.8 


6.9 


3.8 


10.8 




4.5 


4.3 


11.3 


10.5 


3.6 




23420* 


0.126 




0.069 


0.844 


0.057 


0.036 


0.115 




0.045 


0.035 


0.111 


0.079 








12.6 


3.3 


7.5 


81.8 


4.9 


3.3 


10.5 




4.4 


3.5 


12.1 


9.2 






23498 


O.IOO* 






























6.1 x 





























fVariable weather conditions during the exposure of IRAS 16487 in the H band make flux calibration uncertain. 
"(Variable weather conditions during the exposure of IRAS 23420 in the K band make flux calibration uncertain, 
(x) Indicates line detected only at 2cr limit, not above 3<r as is the standard. 
Columns: (1) Name of object. (2) — (15) list the fluxes and equivalent widths for the major lines observed where the top 
value is observed line flux, F (10~ 14 erg s _1 cm" 2 ) and bottom value is observed equivalent width, EW (A). A value in 
italics indicates the lines falls in a waveband where atmospheric transmission falls below 80%, ... indicates a non-detection 
and a blank column indicates the line is at a wavelength rcdshifted out of the observed band. 



IRAS 20210, which is closest to the galactic plane with a 
galactic latitude of -14° and Ay of 0.50 mag. We would ex- 
pect values for Ay using only optical lines to be less than 
those derived from IR lines, since optical li nes probe to a 
smaller optical depth than those in the NIR (Puxley 1991). 



Potential problems with the use of this method exist. 
It is strictly only valid if the dust in these galaxies, known 
to be an important component in ULIRGs, acts as a dis- 
crete absorbing screen in front of the emitting gas. If mixed 
with the gas, dust scattering may have an important effect, 
preferentially enhancing H/3 emission. Absorption features 
from the underlying stellar continuum may have affected 
emission lines for some objects. Observing effects such as 
wavelength dependent slit losses, seeing, guiding effects and 
centring errors may affect the NIR line fluxes, especially im- 
portant here since most of the objects are extended. Values 



for Ha and H/3 fluxes come from observa tions u sing a wider 
slit than used here i.e. Veilleux et al. (1999a) used a 2.0 
arcsec slit, our observations were made with a 1.22 arcsec 
slit width. This could lead to underestimates of NIR line 
fluxes and hence Ay, in cases where the emission line region 
is extended over many arcsecs. Finally, differences in the 
dust coverage could affect these values, potential 'clumpy' 
dust coverage may result in NIR lines coming from a region 
not seen in the optical. Even given all these uncertainties, 
it is expected that general trends in line ratios will still be 
revealed. 



4.3 HI recombination lines 
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Table 4. Values found for extinction using various line ratios in the optical 
and infrared. Note that IRAS 00150 and IRAS 23420 have no optical obser- 
vations available in the literature and hence have not been included here. 



Ay- (mags) 



Galaxy 
(IRAS) 


Ha 
W0 


Paa 
Hp 




Paa 
Ha 


Br'y 
~TTa 


Galactic 1 " (A K ) 


Reft 


16487 


1.49 


3.01 


1.33 


3.78 


1.24 


0.06(0.01) 


(1) 


17179 


3.72 


4.81 


3.81 


5.36 


3.88 


0.09(0.01) 


(1) 


20210 


2.54 


3.19 


2.26 


3.53 


2.08 


0.50(0.06) 


(2) 


23365 


2.79 


4.99 


3.66 


6.14 


4.09 


0.34(0.04) 


(3) 



tTaken from Schlegel et al.( |199&j ) 
fSource from liter ature u sed for optical line ratio . 
References: (1) Veilleux et al. ( 1999a). (21 Y oung et al. (199f), (3) Veilleux 

et al. (1999b|). 



4-3.1 Paa emission 

Several Paschen and Brackett series recombination lines are 
evident in these spectra. The brightest line for each galaxy 
is Paa, at a rest wavelength of 1.876/xm. In every observa- 
tion Paa is narrow, having deconvolved FWHM ~ 300 - 500 
kms~\ except for IRAS 20210 which shows FWHM ~ 700 
kms -1 . Notably, in all but IRAS 00150, a blue asymmetric 
wing at Paa is seen to some degree. This phenomen on ha s 
been noted in ULIRGs previously e.g. Veilleux et al. (1997) 



who find that many of their objects display blue asymmet- 
ric wings at Paa, including IRAS 17179 and IRAS 23498. 
This could be faint broad emission from high velocity gas 
(FWHM ~ 2500 km s" 1 ). However in most cases, with the 
possible exception of IRAS 17179 and IRAS 20210, other ex- 
planations are possible. One possibility is faint H2 emission 
from t he 7-5Q(3)1.8721^m and 6-4Q(5) 1.8665 um transi- 
tions ( |Veilleux et al. 1997| [Veilleux et al. 1999b ). Another 
possibility is the presence of a pair of Hel em ission lines at 
1.8686mhi and 1.8697^m flMurphy et al. 1999| ). At the spec- 
tral resolution of our data we cannot distinguish between 
the possibilities. 



4-3.2 Broadened emission 

The only broad emission observed is Pa/3 with FWHM 
~ 5000 km s~\ seen in the spectrum of IRAS 23498. Previ- 
ous observations of this Seyfert 2 consistently show broad- 
ened emission (FWH M ~ 3000 km s" 1 ) in the Paa line 
([Veilleux et al. 1997). 



4-3.3 L,Br~t VS Lf 



correlation 



A significant correlation between \jbt-, and Lfir has been 
found for LIRGs, the less luminous counterparts of ULIRGs. 
ULIRGs however, do not follow the same rela tionship and 
appear to be underluminous in Br7 emission (Goldader et 



al. 1995). We have added the data for our sample to the 
previously published data to see if they fall upon the same 
relationship. 

For consistency we have calculated Lfi r (= L(40- 
120/im)) using the formula fro m Hel ou et al. (1985). Note 
that, following Goldader et al. (1995), we have taken H D = 
75 km s _1 Mpc -1 for this calculation. Where available, we 
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Figure 3. Dereddened Br7 luminosity vs far infrared luminos- 
ity. Filled circles represent ULIRGs from this sample, filled stars 
represent ULIRGs from this sample uncorrected for extinction. 
Empty squares and circles represent LIR Gs and ULIRGs respec- 
tively from the sample of Goldader et al. (1995). The solid line is 
the best fit to the LIRG sample, also from Goldader et al. (1995). 



have used E(B-V) from the literature, calculated using opti- 
cal line ratios (see Table ^ for more details). Where E(B-V) 
was not available the data were not corrected for extinction. 

Figure ^ shows Ls r7 vs Lfir for both datasets, where 
the solid line is the best fit to the LIRGs made by Goldader 
et al.. They conclude that the ULIRGs from their sample 
fall below the correlation indicating that they are underlu- 
minous in Br7. The two lower luminosity galaxies from our 
sample, classified in this work as ULIRGs, would actually 
be classified in the Goldader scheme as LIRGs, due to the 
different value of H used. Our higher luminosity objects 
do not seem to follow the trend of underluminosity in Br7 
emission with respect to the correlation marked. In fact, our 
data points may follow the correlation if the large scatter 
seen for the LIRGs is taken into consideration. 

The importance of aperture effects must also be con- 
sidered. On average, the objects from the Goldader sample 
have been observed with a 3.0 arcsec wide aperture. The slit 
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width used in our observations was 1.22 arcsec. This may 
result in our observations giving a lower line flux, and hence 
luminosity, by a factor of approximately 2, with respect to 
those of Goldader et al.. 

The depression of the Br7 luminosity in ULIRGs, when 
compared with LIRGs, is explained by Goldader and col- 
leagues by a simple model where ULIRGs contain very com- 
pact nuclei which are optically thick at 2/^m. They reason 
that the Br7 is therefore highly obscured, resulting in a poor 
correlation between Br7 and FIR luminosity. 

An alternative explanation would be to question the 
correlation found for LIRGs. When considered as a group, 
the relationship found can be seen to be very dependent on 
the two highest luminosity points. If these are disregarded, 
then the scatter of the remaining points would change the 
fitting of any trend dramatically. When the additional higher 
luminosity data points provided by the ULIRGs are added 
to the plot then a correlation with a lower gradient could 
be fitted to both sets of data. Indeed, there is such a large 
scatter in the data that it is possible that there is no cor- 
relation and instead simply two 'regions', representing the 
range of properties found for each group. They find support 
for their conclusions in other work in volving LIRGs, where 



similar rela tionships have been found (DePoy 1987; Veilleux 
et al. 1995). The distinction between LIRGs and ULIRGs 



is arbitrary, and claims that they are fundamentally differ- 
ent would need to be supported by observations at other 
wavelengths. 



4-3.4 Second nuclei 

Four of the galaxies observed have multiple nuclei. Of these, 
two were observed with the slit orientation covering both 
nuclei. This was the case for IRAS 00150 and IRAS 20210. 
The spectra of the second nuclei are plotted with the pri- 
mary nuclei, at the same scaling, in Figure |l|. As can be 
seen, the spectra of the secondary nuclei seen are practi- 
cally featureless. This is consistent with preliminary findings 
from another recent NIR survey of 33 ULIRGs (including 
I RAS 16487 and IRAS 23365) conducted by Murphy et al. 
( [l999| ). 

The only feature we see for the secondary nucleus of 
IRAS 00150 is weak Paa emission. The second nucleus of 
IRAS 20210 shows no emission lines, but does show CO ab- 
sorption, indicative of an old stellar population. It is useful 
to make an estimate of the contribution of the secondary 
nucleus to the overall FIR luminosity of the system as de- 
tected by IRAS. For IRAS 00150, assum ing the relationship 
between Lsr-y and Ijfir holds (§ 4.3.3) and using the ra- 
tio of Paa in the two nuclei, we estimate the minor nucleus 
contributes 20 per cent of the FIR emission. For IRAS 20210 
we used the ratio of Br7 emission. An upper limit of Fsr-y 
— 3.81 xW~ 16 erg s~ 1 cm~ 2 was found in the secondary nu- 



cleus, indicating a contribution of ~ 6 per cent. Therefore, 
in both cases, the contribution of the secondary nucleus is 
minor. We note the the similarity of the spectra of the two 
nuclei seen for IRAS 00150, here the minor component could 
be a scaled version of the major. This is not the case for 
IRAS 20210, where the line to continuum ratio seen in the 
major component is not mirrored in the second nucleus. 



4.4 Excitation of H2 lines 

4-4-1 Thermal vs non-thermal emission 

Emission from electric quadrupole transitions within the 
ground state of molecular hydrogen is often bright in in- 
frared observations of starburst galaxies and AGN. There 
are two types of excitation mechanisms producing these 
lines: collisions due to shocks (thermal process) and radia- 
tive decay from an excited state produced by UV photons 
or X-rays, resulting in fluorescence (non-thermal process). In 
general both mechanisms, in varying proportions, contribute 
to the emission. Two galactic cases of particular interest are 
(1) supernova remnants (SNR) which exhibit purely ther- 
mal shock excited emission, and (2) planetary nebulae which 
exhibit pure fluorescent emission. In starburst galaxies OB 
stars produce UV photons, also SNR and winds from mas- 
sive stars drive large scal e outflows which produce shocks 
(e.g. Heckman et al. (199C). Often the more dominant mech- 
anism is uncertain. In AGN, the stronger H2 emitters, the 
emission is believed to be associated with the active nucleus 
which emits a strong UV - X-ray continuum and also drives 
outflows. 

Enhanced H2 emission in mergers may arise from shocks 
produced by cl oud-c loud collisions, as originally proposed by 
Harwit et al. (1987). However the small filling factor of the 
dense molecular clouds presents problems for this model. 
A solutio n to t his problem may be that proposed by Jog & 
Solomon (1992), who describe a scenario in which HI clouds, 
associated with the interacting galaxies and with large filling 
factors, collide. This leads to the formation of hot, ionised, 
high pressure gas which causes a radiative shock compression 
of the denser H2 clouds. 

Since the non-thermal mechanism is more efficient than 
the thermal in the population of v — 2 and higher vibrational 
levels, in principle line ratios can be used to discriminate 
between these mechanisms. 

A minimum of two H2 lines were detected for each 
galaxy, increasing to nine for the brightest galaxy observed, 
IRAS 20210. The potential for deconvolution of this emis- 
sion into thermal and non-thermally excited components is 
small for the galaxies which show only a few H2 lines and so 
only generalised conclusions can be drawn. For the brighter 
galaxies: IRAS 20210, IRAS 23365 and IRAS 23420 we use 
the population and diagnostic diagrams discussed later in 
this section. 



4-4-2 Line ratios: a direct comparison 

First we draw some general conclusions concerning thermal 
vs non-thermal excitation for each galaxy by the direct com- 
parison of observed line ratios with those from models. Ta- 
ble ^ lists the measured ratios compared with typical shock 
and low density fluor escent predictions taken from Black and 
van Dishoeck (1987), plus a mixed model assuming equal 



contributions to 1-0S(1) flux from shock and UV compo- 
nents. 

IRAS 00150, IRAS 16487 and IRAS 17179 have only two 
H2 line detections. This limits our analysis substantially, 
but we can conclude that the ratio is indicative of a shock 
origin over fluorescence. We note that in IRAS 00150 and 
IRAS 17179 the ratio is ~10 per cent less than that expected 
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for shock excitation. However this may be attributed to er- 
rors in the marginal detection of the lines above the noise 
in the spectra. 

Line ratios for IRAS 20210 strongly indicate a purely 
shock origin for the emission with good agreement for all 
but one of the ratios considered. The one discrepancy is the 
2-lS(l) line which was only detected above the 2a level. 

The ratios measured for IRAS 23365 are fitted best by 
the predictions produced by the mixed model. This predicts 
similar fluxes for 1-0S(0) and 2-lS(l), indicating approxi- 
mately equal contributions from shock and fluorescent exci- 
tation processes. 

Finally, considering the ratios seen for IRAS 23420 the 
best model is unclear. The 1-0S(2) line ratio falls between 
predicted values for both the shock and mixed models. The 
low ratio value seen for 1-0S(0) however, the stronger of the 
two lines, favours the shock model over the mixed model. 

As is clear from the above, it is not possible to at- 
tribute the emission uniquely for half of the sample observed. 
Throughout, we are constrained by the lines observed which 
mostly come from the v=\ levels. The population of these 
levels is dominated by thermal excitation which falls off 
rapidly in other bands. Hence to deconvolve this emission 
successfully, detection of lines from the v=2 levels is neces- 
sary. This becomes evident when comparing the predicted 
ratios in Table The best lines to use are those having very 
different ratios for shock and UV models. The 1-2S(1) line, 
which has been detected in every galaxy in this sample, does 
not fulfill this criteria, having only a ~10 per cent difference 
between the models. 2-lS(3) and 2-lS(l) having ~30 per 
cent and 50 per cent change in predicted values are bet- 
ter. Unfortunately they are the weaker lines and are often 
undetected. 



4-4-3 Population and diagnostic diagrams 

Various methods exist to discriminate further between ther- 
mal and non-thermal excitation mechanisms. We use two 
of these methods for the brighter objects in our sample 
(IRAS 20210, IRAS 23365 and IRAS 23420) in order to bet- 
ter constrain the excitation mechanisms. Level population 
diagrams are a useful tool, employing similar theory to the 
previous method, but utilising more emission lines. We also 
employ diagn ostic diagrams with H2 line ratios, used by 
Mouri (1994) to study different thermal excitation mech- 



anisms using characteristic temperatures. 

Le vel population diagrams for TRAS 20210, TRAS 23365 



for both v—1 and 2 levels, as a straight line fit using only 
v=\ lines can be obtained for a galaxy even when a sig- 
nificant amount of fluor escent emission is present, as noted 
by Tanaka et al. ( [l989| ) for NGC7027 and Ml-78, which 
exhibit 10 per cent and 20 per cent non-thermal emission 
respectively. 



Following the method used by Davies et al. (200C), a 
combination of thermal and non-thermal models have been 
fitted. Free parameters were the relative fraction of non- 
thermal emission, fuv , the temperature of the thermal emis- 
sion, T, and the absolute scaling. The thermal component 
was represented by a Boltzmann distribution at temperature 
T given by 

N u -E u 
g u kl 

where E u is the excitation energy of the upper level. The 
non-thermal model used was model 14 from Black and van 



Dishoeck (1987). The best fit was found using a chi-square 



minimalisation technique. 

On scrutiny of Figure ^ it is immediately obvious that 
none of the galaxies exhibits purely thermal emission. The 
data for IRAS 20210, with 9 derived points, indicates a 
mostly thermal origin with a small fraction of fluorescence 
also present. Analysis indicates 10 per cent of the 1-0S(0) 
flux, 40 per cent of the total H2 flux, comes from non-thermal 
emission, with a temperature of 1870K. This model fits well 
for v=\ lines but fits less well for the v=2 lines. Additi onally , 
when placed on the diagnostic diagrams of Mouri (1994), 
IRAS 20210 is in the region occupied by the shock model 
and supernova remnants - objects known to show purely 
thermal emission. 

The population diagram for IRAS 23365 is quite differ- 
ent, a single straight line fit is very unlikely. The lines seen 
for v = 1, J = 6 and 7 and v = 2, J = 3 are weak detections 
and in one case (J=7) only an upper limit, which was not in- 
cluded in the fit. We find for this object a higher proportion 
of non-thermal contribution, 30 per cent of the 1-0S(1) flux, 
70 per cent of the total H2 flux, with a temperature of 2800K. 
This model fits well except for the 1-0S(4) line which falls 
at a wavelength which shows a drop in transmission for this 
object . When plotted on the diagnostic diagrams of Mouri 
(1994), this galaxy falls in the region occupied by NGC 6240. 
NGC 6240 is a galaxy-merger starburst/LINER system with 
much disputed emission mechanisms. Values for the amount 
of non-ther mal H2 emission exhibited by this object range 



and IR KS 23420 are shown in Figure |, derived using the Mouri ( |1994| ) favours shock heating for the thermal 



from none ( 5ugai et al. 1997 ) to 70 per cent (Tanaka et al 



local thermal equilibrium ortho/para ratio of 3. The level 
population is the column density, N„, divided by the statis- 
tical weight, g u , and is calculated using the relationship: 

N u I 



using the measured flux, I, transition probability, A u i, and 
frequency, v. The statistical weight, g u , is given by the prod- 
uct of the rotational and spin degeneracies. The values of 
level population are scaled to that of the v=l, J=3 level. 

Simplistically, if the points on such a diagram can be fit- 
ted by a straight line then the emission mechanism is purely 
thermal, otherwise there is additional emission due to non- 
thermal mechanisms. However data needs to be available 



component of this galaxy. Based on the similar position of 
these galaxies in the diagnostic diagrams we suggest that a 
similar interpretation might apply to IRAS 23365. 

Finally we consider the population diagram for 
IRAS 23420. Unfortunately we detect only lines from the 
v — 1 level for this galaxy which makes the interpretation 
inconclusive. Some of these detections are also only weak, 
leading to uncertainties. Our model finds 30 per cent of the 
1-0S(1) flux, 70 per cent of the total H2 flux, having a non- 
thermal origin, with a temperature of 2450K, agreeing with 
the conclusions drawn previously using the line ratios. 
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Table 5. H2 line ratios compared with predicted ratios. 

Line Wavelength Relative flux of IRAS galaxies* 



Models* 







00150 


16487 


17179 


20210 


23365 


23420 


Shock 


UV 


Mixed 


1-0S(2) 


2.0338 


0.25 


0.32 


0.24 


0.34 


0.42 


0.41 


0.37 


0.50 


0.44 


2-lS(3) 


2.0735 








0.07 






0.08 


0.35 


0.22 


1-0S(1) 


2.1218 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1-0S(0) 


2.2235 








0.25 


0.29 


0.15 


0.21 


0.46 


0.34 


2-lS(l) 


2.2477 








0.13 


0.29 




0.08 


0.56 


0.32 



fNormalised to the flux of 1-0S(1) emission line, 
f Respectively model S2 and model 14 taken from Black and van Dishocck (1987). The mixed 
model assumes 50% of 1-0S(1) emission is thermal and 50% is non-thermal. 



4.5 



Correlation between H2 v = 1-0S(1) and 
[OIU6300 



Mouri et al. (1989) find a linear correlation between the H2 
f=l-0S(l)/Br7 and [OI]A6300/Ha emission line ratios for a 
sample of galaxies containing an AGN or a starburst nucleus. 
Their data show both the H 2 1-0S(1) (hereafter S(l)) and 
[OI] lines are enhanced in AGN. Using optical data, where 
available, from the literature we have added the galaxies 
we have observed to those in the Mouri sample, see Fig- 
ure p. Squares represent starburst nuclei, circles represent 
AGN and triangles represent LINERS. The circled points are 
the four galaxies from our sample: IRAS 16487, IRAS 17179, 



IRAS 2 C210 and IRA3 23305. Classiiicaliuiis fen bfaese galax- 
ies, takl cn from the literature, 



LINER respectively. The positions of these galaxies on this 
plot are as expected; they lie amongst the AGN with en- 
hanced S(l) and [O I] em ission. 

Mouri et al. dl98§) have only one LINER included 



in their sample, NGC 6240 (mentioned earlier) , which they 
note has a remarkable S(l) line width with FWHM ~ 650 
km s _1 . When NGC 6240 is plotted it falls at the extreme 
upper right, away from the other galaxies. The authors inter- 
preted this as a special case in which the emission is excited 
by shocks directly, due to the merger of two galaxies, in 
addition to the effects of a central engine and/or starburst 
activities. This theory has been supp orted by recent work 
on NGC 6240 by Ohyama et al. fl200C|) , who have found ev- 
idence for young superbubbles in the central region. They 
propose that the interaction of this expanding shell with the 
extragalactic molecular gas produces shock-excited intense 
H2 emission. 

The LINER galaxies we have added however, 
IRAS 16487 and IRAS 23365, show smaller S(l) line widths 
(respectively FWHM ~ 400 km s" 1 and ~ 250 km s" 1 ), 
consistent with their position amongst the AGN. 



4.6 [Fe II] emission 

The most prominent [Fell] transition lines seen in the NIR 
are a 6 D 9 / 2 - a 4 D 7 / 2 at 1.257/im and a 4 F g / 2 - a 4 D 7 / 2 at 
1.644/im. Due to the range of redshifts of the objects ob- 
served, the former line was not covered by our data. The 
wavelength of the line at 1.644/wn meant that it could po- 
tentially be observed in every galaxy. Unfortunately, due 
the range of redshifts observed, this line falls in the region 
where atmospheric absorption is high for some objects. This 
is the case for IRAS 00150, IRAS 17179 and IRAS 23420. De- 



spite this the [FeII]1.644^tm line is detected in all the galax- 
ies observed except IRAS 00150, where a slight increase in 
flux can be seen at the appropriate wavelength. Another, 
lower flux, [Fell] emission line at 1.9670^im falls in the wave- 
band observed. This line may be detected in IRAS 17179 and 
IRAS 20210. An increase in flux at the appropriate wave- 
length can be seen for IRAS 23365 and IRAS 23420, however 
there is no detection above the 2 a limit. At this resolution, 
this line would be blended with [SiVI] emission. Such an 
apparently low incidence of the [Fell] 1.9670/im line in this 
sample is perhaps slightly surprising when compared with 
a recent survey of ULIRGs where detections of this line are 



claim ed for 81 per cent of objects observed (Murphy et al 



200C ). However the typical [FeII]/Paa ratio quoted is around 
1 per cent for objects with weaker [Fell] emission, rendering 
the line too faint to be se en for the objects we observe. We 
note Murphy et al. ( ^OOOj ) do not detect [Fell] 1.9670/wn for 
IRAS 16487 or IRAS 23365. 

The excitation mechanism of [Fell] emission is still un- 
der debate. Like [OI], [Fell] emission is excited by elec- 
tron collisions in partially ionised regions (PIRs). Mouri et 
al. (2000) attribute extensive PIRs in starburst galaxies to 



shock heating from supernovae, whereas in Seyfert galax- 
ies they propose these regions arise due to photoionisation 
by nuclear radiation. Alonso-Herrero et al. (1997) demon- 
strate a similar idea, favouring a progression of the propor- 
tion of shock excitation present from HII regions to star- 
burst galaxies, through to Seyferts and supernovae. This is 
demonstrated in the diagnostic diagram they employ where 
the line ratios of [FeII]1.644/Br7 vs [OI]6300/Hq produce 
a separation of the galaxy types, indicating the proposed 
progression. 

When the objects observe d here are added to the fig- 
ure from Alonso-Herrero et al. (1997) they fall amongst the 
other galaxies, in the area marked by the authors as the 
overlap between Seyferts and composite galaxies. From this 
we would infer a mixture of shock heating and photoion- 
isation as the mechanism for the production of [Fell] in 
these galaxies, with the latter process appearing to domi- 
nate. Overall, based on these ratios the emission seems to 
be typically Seyfert-like, as expected for these galaxies. 



4.7 [SiVI] emission 

Detection of [SiVI] emission indicates the presence of high 
energy ionising photons. This high excitation line has an 
ionisation potential of 167 eV and cannot be formed by con- 
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Figure 5. A plot of H 2 v = 1-0S(1)/Br7 and [OI] A6300/Hq. 
Symbols are as follows: squares represent starburst nuclei, circles 
represent AGN and triangles represent LINER S. Cir cled points 
are from this work, all others from Mouri et al. (h98£|). Solid line 
represents best fit to starburst nuclei data. 



tinuum emission from young stars. It is a feature only associ- 
ated with active galaxies, suggesting that the mechanis m for 



its production i s photoionisation by an active nucleus ( Mar- 



coni et al. 19941 ) . The [SiVI] line is detected for IRAS 17179 
and the main nucleus of IRAS 20210, having FWHM ~ 1200 
km s _1 in both detections. This detection confirms the sta- 
tus of these two galaxies as Seyferts. [SiVI] is not detected 
in the spectra of the LINER type galaxies. The absence of 
this line in the spectra of the galaxies of undefined type: 
IRAS 00150 and IRAS 23420, does not support the presence 
of a Seyfert nucleus, however the presence of a highly ob- 
s cured AGN cannot be ruled out using this criteria alone 
( |Marconi et al. 1994| ). 



4.8 CO absorption 

CO absorption arises from the presence of older stars, red 
giants and supergiants. There are three main bandheads 
potentially visible: 12 CO(2,0) at 2.294/im, 12 CO(3,l) at 
2.323^m and 13 CO(2,l) at 2.345A«n. In our observations, 
due to the redshift range of the objects observed, detec- 
tion of these absorption features was potentially only possi- 
ble in two objects: IRAS 20210 and IRAS 23365. Both ob- 
jects show absorption, possibly visible in both nuclei for 
IRAS 20210. Figure ^ shows an expanded view of the spec- 
tra of the objects in the relevant waveband, useful to show 
the CO absorption not obvious with the scaling in Figure |l| 
The equivalent width of the absorption in the main nucleus 
of IRAS 20210 is smaller than that of IRAS 23365. This is 
as expected, since emission from heated dust begins to be- 
come significant around this wavelength, rising into the mid- 
infrared. In IRAS 20210, a Sy2, the stronger emission due 
to hot dust effectively dilutes the measured CO absorption 
bands. In IRAS 23365 the hot dust continuum is less, allow- 
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Figure 6. Section of K band spectra for IRAS 20210 (both nuclei) 
and IRAS 23365 showing CO absorption in more detail. Note: 
objects are all plotted 'as is' and have not been scaled for viewing 
purposes. 



ing CO absorption to become more evident. In the second 
nucleus of IRAS 20210, labelled in Figure | as IRAS 20210b, 
the emission is characteristic of an older stellar population: 
showing CO absorption and a lack of recombination emission 
lines. Additional evidence for the presence of a stronger 'di- 
luting' hot dust continuum comes from the relatively redder 
spectrum of IRAS 20210 when compared with IRAS 23365. 



5 CONCLUSIONS 

We have presented NIR observations of a sample of mainly 
interacting ULIRGs, comprising of H and K band spectra. 
From the analysis of these data our main aims were to inves- 
tigate the power source of these extremely luminous objects, 
the potential excitation mechanisms of the strong molecular 
hydrogen emission which is often seen in such objects, along 
with an investigation of the other characteristics such as 
extinction and multiple nuclei. Our results are summarised 
below. 

Of the seven galaxies observed, five have previous 
classifications and the remaining two, IRAS 00150 and 
IRAS 23420, have no published optical data other than that 
used for redshift determination. Broadened HI recombina- 
tion emission was detected in one object, IRAS 23498, with 
FWHM ~ 5000 km s _1 . This finding is consistent with pre- 
vious results for this object which is classified as a Sy2. The 
high excitation line [SiVI] is detected in two galaxies in our 
sample, both are classified as Sy2s. We have investigated 
several proposed correlations, which as well as exhibiting a 
specific trend between line ratios, also help to distinguish 
AGN from starburst galaxies. In most cases the galaxies 
added to the existing data exhibited the expected trend . 
We question a correlation found by Goldader et al. ( 1995j ) 
between Ls r7 vs Lfir for LIRGs. This is due to the small 
size of the sample on which the relationship is based. Data 
points added for the ULIRGs suggest either (1) the need 
for a different form of correlation, or (2) that no significant 
correlation is seen. 



Regarding the classification of IRAS 00150 and 
IRAS 23420, neither spectrum exhibits [SiVI] emission. How- 
ever this does not rule out the presence of highly obscured 
AGN. It would be useful to obtain good quality optical data 
for these objects, in order to classify them. 

Another objective was to provide information on the 
excitation mechanism of the molecular hydrogen emission 
observed in the K band. Due to the variation in the number 
of lines observed, several diagnostic methods were used. 

The overall finding for the galaxies in our sample, with 
the possible exception of IRAS 23365 and IRAS 23420, is the 
predominance of thermal excitation. Non-thermal emission 
is also present in some objects but in varying proportions. 
Comparison of only a few line ratios indicated very little, and 
possibly negligible, contribution from non-thermal emission 
in IRAS 00150, IRAS 16487 and IRAS 17179. More detailed 
analysis was possible for the other objects in the sample. 
IRAS 20210 shows only a small proportion, 10 per cent, of 
non-thermal contribution to the 1-0S(1) emission. In con- 
trast, values of 30 per cent non-thermal emission were found 
for both IRAS 23365 and IRAS 23420. This corresponds to 
~70 per cent of the total H2 emission. 

Since most of the galaxies observed are interacting, one 
might expect thermal excitation through shocks to be the 
dominant mechanism. The obvious exceptions to this are 
IRAS 23365 (which indeed shows a higher proportion of non- 
thermal emission), and also IRAS 17179 and IRAS 23420 
which are classified as non interacting galaxies.We note that 
for IRAS 17179 only a few H2 emission lines were detected, 
so our conclusions are less certain. However, the increased 
fraction of non-thermal emission in IRAS 23420 is consistent 
with it being non-interacting. 

Values for reddening were derived for galaxies with 
available optical data, using NIR emission lines in combi- 
nation with the optical lines. The results depended on the 
combination of lines used, with values using only optical 
lines being consistently smaller than those found using a 
NIR line with an optical line. This is as predicted since the 
infrared lines can be detected from more highly reddened 
regions. 

Of the four galaxies with optical data, IRAS 23365 
showed the highest values of extinction with an average value 
of Ay = 4.7 ± 0.8 mag. In comparison the value found using 
optical lines only was Ay = 2.8 mag. The former value is 



consistent with that derived by Davies et al. (2000b), Ay 



= 4.4 mag, found using a model of 500K dust contributing 
50 per cent of the K band emission using the NIR imaging. 
The values of extinction found for the other galaxies were 
moderate, being in the range Ay = 2.3 - 4.5 mag. 

There is a bias in our sample towards interacting sys- 
tems. For the two for which we have spectra for both nuclei, 
we find either very weak or no emission lines, indicative of 
an evolved population, in contrast to the brighter primary 
nucleus. This trend has also been found for other ULIRGs 



( [Murphy et al. 1999j ; [Murphy et al. 2000[ ). An estimate of the 
contribution of the secondary nucleus to the IR luminosity 
of the system was made and in both cases the contribution 
was found to be minor. Observations of the secondary nu- 
clei of IRAS 16487 and IRAS 23498 would be interesting, to 
see if these characteristics are shared by secondary nuclei of 
interacting ULIRGs in general. 

As mentioned above, all the hydrogen recombina- 
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tion lines detected were narrow, with the exception of 
IRAS 23498. An interesting feature that arises frequently in 
these objects is the blue shoulder of the Pacv line. This fea- 
ture could be explained by the presence of a blue velocity 
component. Another possibility is the presence of blended 
lines, as seen at our spectral resolution, higher resolution 
spectroscopy is required to make further progress. 
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APPENDIX A: DETAILS FOR INDIVIDUAL 
GALAXIES 

A short synopsis is given here for each galaxy observed. It 
will include a mention of previous work and a brief summary 
of the emission lines seen. 



Al IRAS 00150+4937 

Little previous data exists for this double nuclei galaxy 
which means it is unclassified. It is amongst the faintest 
of the objects in the sample with small EW for the lines 
detected. It is the only object in the sample for which no 
blue wing of Pacv was seen. Detection of only two H2 emis- 
sion lines made the determination of excitation mechanism 
difficult. The ratios seen were indicative of shock excitation 
over fluorescence. The spectra for the secondary nucleus of 
this system showed very little, except weak Paa emission. 
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A2 IRAS 16487+5447 



This object was classified by Leech et al. (1994) as isolated 
and undisturbed o n the basis of photographic data. Fur- 



ther NIR imaging (Clements & Baker 1996b; Murphy et al 



199f) has shown this object is interacting, although we only 



present spectra of the primary nucleus. The optical classifi- 
cation of this object is LINER , while the ISO classification 
is starburst (Lutz et al. 1999). The Paa line for this ob- 
ject shows a blue wing. Only two Eh emission lines were 
detected, from these we conclude pure shock excitation is 
the most likely excitation mec hanism , consistent with the 
conclusions of Taniguchi et al. ( |l999| ). 



A3 IRAS 17179+5444 



This S eyfert 2 ( Veilleux et al. 1999a) was classified by Leech 
et al. (1994) a s being a interacting system. This work and 



Davies et al. (2000b) show the secondary nucleus to be a 
foreground star. This galaxy is one of the weaker of the sam- 
ple and its classification as a ULIRG depends on the value of 
H used. Paa shows a blue wing which could conceivably be 
broadened emission. [SiVI] is detected for this galaxy. Only 
two H2 emission lines were detected, from these we conclude 
shock excitation is the most likely excitation mechanism. 



lines were detected to allow some analysis of the potential 
excitation mechanism. Final conclusions attributed perhaps 
as much as ~ 70 per cent of the total H2 emission to non- 
thermal mechanisms. Fairly strong CO absorption bands are 
observed are observed in this object. 



A6 IRAS 23420+2227 

Little previous work and no classification exists for this rel- 
atively faint galaxy which shows no sign of interaction. The 
spectra shown here do not support the possibility that this 
galaxy containing an AGN. The Paa emission seen is nar- 
row, possibly showing a blue wing. Three H2 lines are identi- 
fied which indicate a mainly thermal excitation mechanism. 



A7 IRAS 23498+2423 



This interacting Seyfert 2 (Veilleux et al. 1999a) has the 
highest redshift of our sample. We obtained only a H band 
s pectrum for this obj ect since K band spectra already exist 
(Veilleux et al. 1997). Consistent with this previous work, 
our spectrum shows broadened Pa/3 emission FWHM ~ 5000 
km s _ . 



A4 IRAS 20210+1121 



This interacting Seyfert 2 ( Perez et al. 199"c| ; Vader et al 



1992 ) has been relatively well studied in the past. The spec- 



tra we have obtained for this system are the best of the sam- 
ple with high S/N allowing identification of many emission 
lines, included fainter H2 excitation lines. Data in the opti- 
cal, NIR, radio and X-ray exist for this object. In this work 
this object shows the broadest Paa emission of the sample 
with FWHM ~ 700 km s -1 . This line also shows a blue 
wing which could conceivably be further broadened emis- 
sion. Analysis of the H2 emission points towards a mainly 
thermal origin with a fraction of fluorescence. [SiVI] emission 
is detected from this galaxy consistent with its classification 
as a Sy2. Spectra for the secondary nucleus of this system 
were also obtained for this work. These spectra show no 
emission lines. CO absorption bands were detected for both 
nuclei. 



A5 IRAS 23365+3604 

This single nucleus o bject is generally classified as a LINER 
( Veilleux et al. 1995) but has a lso been labelled as a com- 
posite galaxy ( ]Baan et al. 1998 ). It is one of the objects in 
this sample classified as non-interacting but does show a dis- 
turbed morphology. Previous K band spectra exist for this 
object flGoldader et al. 1995| ; [Veilleux et al. 1999b| ). Unfor- 
tunately removal of atmospheric effects was not totally suc- 
cessful for this object, making the spectrum noisy between 
~ 1.82/im and 1.95^tm (observed wavelength). However this 
region was free of expected emission due to the redshift of 
this object. Of the sample, IRAS 23365 showed some of the 
highest values for E(B-V) calculated using NIR lines. The 
Paa emission measured showed a blue wing seen on some 
many of this objects. Th is line was not broad ened, consis- 
tent with previous work ( Veilleux et al. 1999b ). Enough H2 
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